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Topological chiral phases are ubiquitous in the physics of the Fractional Quantum Hall Effect. 
Non-chiral topological spin liquids are also well known. Here, using the framework of projected 
entangled pair states (PEPS), we construct a family of chiral spin liquids on the square lattice 
which are generalized spin-1/2 Resonating Valence Bond (RVB) states obtained from deformed 
local tensors with d-\- id symmetry. On a cylinder, we construct four topological sectors with even 
or odd number of spinons on the boundary and even or odd number of (Z 2 ) fluxes penetrating 
the cylinder which, we argue, remain orthogonal in the limit of infinite perimeter. The analysis 
of the transfer matrix provides evidence of short-range (long-range) triplet (singlet) correlations as 
for the critical (non-chiral) RVB state. The Entanglement Spectrum exhibits chiral edge modes, 
which we confront to predictions of Conformal Eield Theory, and the corresponding Entanglement 
Hamiltonian is shown to be long ranged. 

PACS numbers: 75.10.Kt,75.10.Jm 


Introduction - The study of topological states of mat¬ 
ter is the focus of research of a multidisciplinary scien¬ 
tific community. On the one hand, they give rise to very 
exotic behavior associated to the topological order [T], 
something which has no analog in standard Condensed 
Matter Physics. On the other, it provides an alterna¬ 
tive and promising approach to fault tolerant quantum 
computing [ 2 ]. Furthermore, experiments both with new 
materials [3] and cold atoms [4] are now at the position 
of engineering some of those states, and thus provide a 
laboratory to verify the theoretical findings. 

Two of the main theoretical challenges in the study of 
topological matter are the identification of simple mod¬ 
els and descriptions which can lend us a complete under¬ 
standing of their exotic properties, as well as the deriva¬ 
tion of parent Hamiltonians, so that we get a hint of 
how to create the states, or under which conditions we 
can find them at low temperature. A possible approach 
to attack both challenges is through the use of tensor 
networks, and in particular of Projected Entangled Pair 
States (PEPS) [SJ [ 6 ]. They provide a simple local de¬ 
scription of many body states from which their physical 
properties can be determined, and they imply the exis¬ 
tence of a parent local Hamiltonian. 

So far, most of the PEPS studies have concentrated 
on non-chiral topologically ordered states, like the toric 
code 12 E], double and Levin-Wen models [7], or the 
short-range spin- 1/2 RVB states [SUTT]. The latter is 
specially appealing as it respects su( 2 ) symmetry and 
is translationally invariant, something which naturally 
appears in many materials [ 12 ]. Furthermore, in the 
Kagome lattice it corresponds to a Z 2 spin liquid mm, 
with very well known anyonic excitations. 

Another class of topological ordered matter corre¬ 
sponds to the Fractional Quantum Hall State (FQHS) 


exhibiting protected chiral edge modes m that are fully 
characterized by a chiral CFT El- For example, the 
edge states of the bosonic Laughlin m, Moore-Read m 
and Read-Rezayi m states can be described by a chi¬ 
ral SU{2)k CFT of the Wess-Zumino-Witten (WZW) 
model El, where k = 1; 2 and 3 respectively. The edge 
modes are protected by the long-range topological or¬ 
der of the bulk which hosts quasiparticles (or anyons) 
with Abelian {k = 1 ) and non-Abelian {k > 1 ) frac¬ 
tional statistics. Recently, an emergent topological chi¬ 
ral spin liquid (CSL), similar to the one proposed by 
Kalmeyer and Laughlin (KL) [181 HI] 5 has been discov¬ 
ered in the Kagome quantum antiferromagnet with bro¬ 
ken time-reversal (TR) symmetry [20] as well as in the 
square lattice [ 21 ], realizing the spin analogue of the 
bosonic Laughlin state. A FQHS was also found in a 
TR invariant Kagome Heisenberg model [22]. Whether 
PEPS can be derived for gapped chiral phases of inter¬ 
acting systems remains a major and subtle issue debated 
in recent papers [23l |24]. Nevertheless, a critical chi¬ 
ral PEPS with topological order has been constructed 
by Gutzwilier projecting two copies of chiral states of 
free fermions [25]. Using the PEPS bulk-edge correspon¬ 
dence, the edge CFT could be identified as SO{2)i but 
with an infinite correlation length. 

In this paper we show that by deforming the tensors 
corresponding to the RVB spin-1/2 state, implementing 
d P id symmetry, one can break the TR symmetry but 
keeping the su( 2 ) and other symmetries of the state, and 
obtain a chiral topological state. We analyze the topo¬ 
logical degeneracies, the entanglement spectrum and the 
entropy via the bulk-boundary correspondence. The chi¬ 
ral edge modes bear some features of CFT although a 
precise identification to simple known CFT was not pos¬ 
sible. We work on a square lattice, and find that, as 
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the RVB state |26l [27] , the family we construct has infi¬ 
nite correlation length and, hence, moves along a critical 
phase of a local Hamiltonian. 

PEPS framework-We consider a two-dimensional ( 2 d) 
bipartite square lattice where a rank-5 tensor of 

dimension is assigned to each site. Here, the index 
5 = 0,1 stands for the local physical spin- 1/2 degrees 
of freedom (of dimension d = 2) and the subscript in¬ 
dices l^u^r^d label the virtual states on the four bonds 
as shown in Fig. [^a). Basically, we can write a general 
spin- 1/2 ansatz wave function on a periodic manifold of 
M sites as l^') = |si> ^2 • • • sm), and 

Csi,s 2 ...sm = Contractwhere all tensors 

share the same bond variables with their 4 neighbors and 
“Contract” means that one sums up over all bond vari¬ 
ables of the 2d tensor network. In order to construct 
a generalized spin-singlet RVB wavefunction, we assume 
that the virtual states belong to the 1/200 spin repre¬ 
sentation of dimension D = 3 (the virtual states | t), \ i) 
and | 0 ) carry Sz = 1 / 2 ,— 1/2 and 0 , respectively) PITO]. 
Since in each configuration of the RVB state, exactly one 
bond singlet | fi) ~ | it) is connected to each lattice site, 
the precise sign structure of the RVB state depends on 
the orientations of the latter : hereafter singlets are ori¬ 
ented from one sublattice (A) to the other sublattice (B) 
as shown in Fig. [^b). After a 180-degrees spin-rotation 
on the B sites (see Fig. [^c)) the RVB state takes the 
form of a translationally invariant PEPS with the same 
tensor on both sublattices. Whenever open boundaries 
are present, one needs to assign the values of the non- 
contracted virtual indices at the boundaries to fully char¬ 
acterize the PEPS wave function. 

Chiral RVB - Our goal is to find a PEPS representation 
of a family of states 4^ with su(2), 90 degree rotation 
and translational symmetry, and that contains the RVB 
state on a square lattice. More precisely, let us define 
by Gi the group transformations of (^ 4 ^ corresponding to 
reflection along the vertical, horizontal, and 45 degrees 
axis. We want to impose that coincides with |4^*). 

Eor PEPS, one can enforce this symmetry at the level of 
the tensor A itself. This constraint, together with the 
additional ones imposed by the singlet character of the 
PEPS fully determine the tensor A in terms of three real 
parameters (see Appendix): 

A = AiRi 0 A2R2 0 '<^Achiral ^5 (1) 

where the (real) Ra and / tensors transform according 
to the Bi and B2 irreducible representations (IRREP) of 
C4V [28l|29], respectively. In other words, the A tensor 
has dp id (meaning dx 2 _y 2 0 i dxy) symmetry. On a fi¬ 
nite L X L torus with an even number of sites, the real 
(imaginary) parts of the PEPS are sums of products of 
even (odd) numbers of R by even (odd) numbers of I ten¬ 
sors. Therefore, the PEPS takes the form, T = 
where the real components and transform accord¬ 
ing to the Ai ( 5 -wave) and A2 (g-wave) IRREP, respec- 



FIG. 1: (a) Local tensor A. (b) Square lattice with (nearest- 
neighbor) singlet bonds oriented from sublattice A to sub¬ 
lattice B. (c) Under a 180-degrees spin rotation on e.g. the 
B sublattice, oriented singlets are transformed into symmet¬ 
ric I tt) 0 I 44) maximally entangled pair states, (d) Transfer 
matrix on a periodic ring connecting the (non-contracted) BK 
virtual indices on the left of the ring to the ones on its right. 


tively, so that and Gf^g) = —14^^)- Eor 

A 2 = Achiral = O 5 one recovers the simple nearest-neighbor 
RVB state [9l |T0] . A 2 and Achiral introduce longer range 
(AB) singlet bonds via “teleportation” as in the non- 
chiral spin liquid of Ref. [TTl Note that we found that 
|T^) 7 ^ 0 only if A 2 7 ^ 0 and Achiral 7 ^ 0 simultaneously. 
Hereafter we fix A 1 = A 2 = 1. 

In order to identify the above PEPS as a chiral spin 
liquid with topological order we have to prove that (i) T 
and T* are different in the thermodynamic limit and (ii) 
there are different topological sectors on a cylinder e.g. 
with even or odd number of spinons on the boundaries or 
with even or odd number of (Z 2 ) fluxes penetrating the 
cylinder, (iii) In addition, we have to fully characterize 
the edge physics using the PEPS bulk-boundary corre¬ 
spondence, so that we can determine the Entanglement 
Spectrum (ES) and the Entanglement Entropy (EE). 

Transfer matrix and topologieal seetors - To do so we 
place the PEPS wavefunctions T and T* on a (horizon¬ 
tal) cylinder of circumference Ny (= 4, 6 , 8 ) and length 
Nh 00 . Their overlap can be written as a bra- 
ket (BK) product (T|T*) = Vieftf^^Vight where f is 
a X transfer matrix (TM) and Veft and Vight 

stand for left and right vectors (of dimension ) defin¬ 
ing the boundary conditions. To construct F one first de¬ 
fines a local rank-4 E tensor by contracting the physical 
index of two superposed bra and ket A tensors, namely 
^LURD = ^lurd^i'u'r'd' where the indices of E com¬ 
bine the indices of the superposed bonds of the top (bra) 
and bottom (ket) A tensors. Then, one builds a (ver¬ 
tical) periodic array of Ny such E tensors, contracting 
over the (vertical) bond indices as shown in Fig. [^d). 
Similarly, one can also construct the ’’regular” transfer 
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matrix T as Elurd = 'Es^iurdi^i'u'r'd')* which 
the normalization (^|^) can be obtained. If the lead¬ 
ing eigenvalue (LE) 7ee of T (normalized by the LE of 
r) is strictly smaller than 1, then ^ and are orthog¬ 
onal in the oo limit. This is what is observed 

in Eigs. I^a) as soon as T acquires an imaginary com¬ 
ponent. Also, for increasing circumference, jee ^ 0 as 
seen in Eig.j^b), implying that T and T* are indeed two 
independent states. 

Topological properties of T (topological sectors, etc...) 
can be obtained from the “regular” TM [30]. It is easy 
to check that T (like T) has a simple block-diagonal 
structure associated to a number of conserved quanti¬ 
ties like the difference — 5'^®^ between the 

2;-components of the total virtual spin m on the bra 
and the ket, namely T = 0r[5'^^]. We found that 
the leading eigenvalues of the TM’s appear in the low¬ 
est sectors, namely = 0 and = ±1/2. 

Secondly, one finds that the parity of the number of vir¬ 
tual |0) states (or spinons at the boundaries) is conserved 
independently in the bra (top) and ket (bottom) layers, 
as for the Kagome RVB Z2 liquid |9l[l0l|32], suggesting 
topological order. The = 0 sector of T is therefore 
split into even-even (ee) and odd-odd (00) sub-sectors, 
namely r[0] = Tge 0 Too, while the = 1/2 sector 

is split into two degenerate even-odd (eo) and odd-even 
(oe) sub-sectors, namely r[l/2] = Tgo 0 ^oe- addi¬ 
tion, one can insert (horizontal) strings of Z2 vison oper¬ 
ators IS HO] in the bra and/or ket layers of T (restricting 
here to Tee)- We found that the resulting TM depend 
only on the parities of the number of Tr-fluxes in each 
layer. We denote by Tt^q (rTTTr) Ihe TM with a vison flux 
in a single (both) layer(s). 

Erom the above analysis we see that, a priori, one 
can construct four wave functions TJ and TJ 

characterized by the parities of the numbers of spinons 
at the boundary (as specified by the subscript) and 
fluxes penetrating the cylinder (as specified by the super¬ 
script). However, we still have to check that these wave 
functions, locally indistinguishable, remain truly distinct 
(and hence orthogonal) in the thermodynamic limit. Eor 
Nh > Ny, the overlaps behave as (T^|Tg) = (feo)^^ 
and with feo = leo!{leelooY^'^ and 

/ttO = 77ro/(7ee77r7r)^/^, where 7ee, 7oo, 7eo, and 7^0 

are the LE of Tge, Too, Tgo, and respectively. 

If the wave function subspace remains degenerate in the 
Ny ^ oc limit, we expect that (i) ^00 and 77^7^ converge 
exactly to 7ee (the largest LE set to 1) while (ii) feo and 
/77O converge to values strietly smaller than 1. In con¬ 
trast, if /t^o ^ 1 (or feo 1)5 the insertion of a flux 
in the cylinder (or a spinon at the boundaries) does not 
generate a new topological sector. The finite size scal¬ 
ings shown in Eig. [^b) - or in the Appendix for other 
parameters - are all compatible with (i). However, ex¬ 
trapolations of either / t^q or feo might be inaccurate in 
some cases - see Eig. [^b) and Appendix for other pa- 



FIG. 2: (a) Leading eigenvalue ±e (normalized by 7ee) of the 
Lee block of the transfer matrix as a function of Achiral (for 
Ai = A 2 = 1) and for At; = 4, 6 and 8 (as shown on plot), (b) 
Finite size scaling of the leading eigenvalues of the diagonal 
blocks (Fee and Foo) and of the off-diagonal blocks (Fee and 
Feo) of the transfer matrix (after proper normalization). The 
diagonal (off-diagonal) blocks have = 0 = 1/2) 

quantum numbers, (c) Finite size scaling of the second leading 
eigenvalues of the Fee, Foo and Feo blocks (after proper 
normalization). 

rameters - due to slow convergence and the limitation to 
small perimeters. It is however plausible that feo < 1 
and /770 < 1 in the thermodynamic limit (at least for 
some parameters) ensuring the existence of four Inden- 
tity (Tg), Spinon (T^), Vison (Tg) and Vison/Spinon 
(Tg) topological sectors (wave functions). 

Asymptotie eorrelations - The behavior of the second 
LE (to be normalized by the leading one) of each 
block of F, provides key informations on the long¬ 
distance correlations associated to the quantum num- 
ber of the block - typically one expects 'iff ^ 
exp (—l/^^jy). The finite size scalings shown in Fig. j^c) 

suggest that both 7ee^ ^ 1 and ^^o h 00 1- We also 

find that the next LE are also going to 1, an evidence 
of gapless even-even and odd-odd sectors of the TM i.e. 
Cee = Coo = c>c. Correlations of spin-singlet operators 
are therefore expected to be algebraic (infinite correla¬ 
tion length). In contrast, we find a gapped even-odd 
sector with an extrapolated value 7oe/Toeloo — 0-85 cor¬ 
responding to a magnetic correlation length Ceo ^ b (in 
units of the lattice spacing) for Achiral = 1- 

Entanglement Speetrum and ehiral edge modes - Us¬ 
ing the standard procedure [6|, we have computed the 
ES, EE and Entanglement Hamiltonian (EH) for infinite 
cylinders (of circumference Ny =4,6 and 8) cut into 
two (semi-infinite) halves. The fixed-point right and left 
vectors, and obtained by applying the TM it¬ 
eratively on (random) boundary vectors belonging to the 
= 0 even (/i = e) and odd (/i = o) sectors are 
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FIG. 3: Re-scaled Entanglement Spectrum vs momentum 
along the edge (for Ai = A 2 = Achiral = !)• (a), (b) and (c) 
show different sectors of the Z-component of the virtual spin 
on the edge. The dashed lines show the chiral edge modes. 
In (a), the momentum of the states of the chiral mode is 
defined mod tt. The multiplicities 1, 1, 2, and 3 of the quasi¬ 
degenerate states in the boxes - expected for a GET - are 
shown on the plot. 

viewed as operators (of trace 1) acting on the virtual 
indices at the edge. From the right and left operators 
aL = + rr) and an = \{V^^ + one obtains 

the EH He = He is block- 

diagonal and the blocks are labelled by the momentum K 
along the edge and by the modulus of the z-component 
of the virtual (edge) spin |5'^| = [^. We 

have computed the ES on infinite cylinders as a function 
of K. At low (quasi)energy, linearly dispersing chiral (i.e. 
with a definite sign of their velocity) modes with almost 
equally spaced levels are seen for all values of Ny (see 
Appendix). Fig. shows the Ny = S spectra which have 
been shifted and rescaled (by the same amount) to set the 
average level spacing of the edge states to A = 1. This 
makes the resemblance with the spectrum of a chiral CFT 
very striking, where each “tower of states” corresponds 
to a WZW primary field and its descendants. Since the 
ES inherits from the spin-singlet character of the PEPS 
SU(2) Kramer degeneracies, the SU{2)i CFT of the KL 
state is the most natural candidate. However, the SU (2)i 
CFT has only two primary fields - the Identity and the 
Spinon fields - in disagreement with the plausible ex¬ 
istence of four topological sectors for the d ^ id RVB. 
In addition, our numerical estimation of the conformal 
weight hi /2 of the Spinon sector (from the ratio of the 
lowest \Sz\ = 1/2 energy level to the \Sz\ = 1 one) gives 
hi /2 — 0.93 for Ny = 6 and hi /2 — 0.84 for Ny = 8, 
significantly above the value hi /2 = 1/4 for SU{2)i. 

Entanglement Entropy and Hamiltonian - We have 
computed the Von Neumann EE from the ES and re¬ 
sults for the two topological sectors (with and without 


EIG. 4: (a) Entanglement entropy S vs cylinder perimeter in 
the two topological sectors with (open symbols) and without 
(full symbols) Z 2 flux (for Ai = A 2 = Achiral = !)• Linear fits 
(area law) S = aNy + b are shown, (b) Weights vs N of the 
Entanglement Hamiltonian decomposed in terms of A-body 
contributions, for Ny = 6 and Ny = S. 

Z2 flux) are shown in Fig. |^a). In the Identity (integer 
spin) sector and no flux one can fit the EE according to 
the area law with a negative intersect b ^ —0.596. This 
value does not agree with the expected — In 2/2 value for 
SU{ 2)1 but is close to —In2 ^ 0.69 also found in the 
critical RVB spin liquid m- To investigate the range 
of the EH we have decomposed it in terms of separate 
A-body terms whose weights can be computed. As can 
be shown in Fig.|^b), the BH is long-range and retains a 
large weight on extended operators involving all the sites 
of the ring i.e. with A = Ny. 

Diseussion and outlook - It is well known in Particle 
Physics and lattice gauge theories (Nielsen and Ninomiya 
no-go theorem [33]) that a chiral field theory cannot be 
discretized on a lattice. This has brought some doubt 
in the community about the possibility of constructing 
a PEPS (of finite bond dimension D) that could be the 
ground state of a local Hamiltonian within a topologi¬ 
cal chiral phase. Together with Ref. |25l the present work 
brings new perspectives : we have constructed a (two pa¬ 
rameter) family of PEPS which are CSL with topological 
order and chiral edge modes and share some properties 
with the KL spin liquid state. Nevertheless, a precise 
identification of the PEPS was not possible and some 
differences with the KL state may arise : (i) First, once 
the time reversal is broken, the KL state is expected to 
have two-fold topological degeneracy on a torus. A pri¬ 
ori, our CSL bears four topological sectors. The scal¬ 
ing of the relevant leading eigenvectors of the TM does 
not allow to safely conclude that two sectors only will 
survive in the thermodynamic limit, (ii) Our numerical 
simulations also point towards critical (singlet) correla¬ 
tions. This property may be closely related to the bi- 
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partiteness of the lattice and the existence of an effective 
field theory formulated in terms of height/gauge degrees 
of freedom [34]. Similar construction on non-bipartite 
lattice like the Kagome lattice (requiring larger bond di¬ 
mension D > 3) is left for future studies, (iii) Lastly, a 
precise assignment of the edge physics to a simple CFT 
was not possible, also calling for further studies. Devi¬ 
ations from a simple SU{2)i CFT (expected for the KL 
state) may appear due to the critical property of the bulk. 
Finally, we note that PEPS construction of non-Abelian 
CSL with edge modes described by SU{2)k CFT with 
k> 1 would also be of great interest. 
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Appendices 


DERIVATION OF THE PEPS ANSATZ. 


Here we wish to construct a modified RVB complex singlet wavefunction whose real and imaginary components 
transform according to different irreducible representations (IRREP) IZi and IZ 2 of the (^4^ point group of the square 
lattice [28]. More precisely, (i) the wavefunction should have the form |^) = |^7^l) +i|^7^2): where |^7^l) and |^7^2) 
are assumed to be real (and non-zero), and (ii) IZi and IZ 2 have to be such that, under any reflection symmetry of 
Cavi I^) transforms into |^*), its time-reversed state (possibly up to a sign), which is indeed a necessary condition for 
a chiral spin liquid. With PEPS, one can achieve exactly this goal by enforcing symmetries at the level of the tensor 
itself. Let us write A = + iJ, where R and I are two real tensors. We impose that R and / transform, under point 

group operations, according to the Ti = Bi {dj. 2 _y 2 orbital symmetry) and T 2 = B 2 {dxy orbital symmetry) IRREPs, 
respectively. In that case, it is easy to see (see main text) that the real and imaginary parts of the wavefunction (on 
a L X L torus) transform according to the IZi = Ai (s-wave) and IZ 2 = A 2 (g-wave) IRREP of Cav 
Hence, we require that the R and / components of the A tensor follow. 


Tys _ JDS _ JDS _ _ JDS _ _ JDS 

-^lurd -^Idru ^ruld ^drul -^uldri 

JS _ _ JS _ _ JS _ JS _ JS 

^lurd -^Idru ^ruld ^drul -^uldr' 


( 2 ) 

(3) 


Here l^u^r and d label the bond variables (0,1, and 2 for the | f), | i) and |0) virtual states, respectively) clockwise 
around the site, starting from its left bond, and s = 0,1 is the physical index. To solve these equations one has to go 
through all possible combinations of the bond indices /, r and d. One obtains. 


In fact, spin SU(2)-invariance implies that. 


^2s22 

\ ss 

— : 

(4) 

^222s 

\ ss' 

— '^1 5 

(5) 

^22s2 

\ ss' 

— 5 

(6) 

■^s222 

= -Af . 

(7) 

: 

— ^ss' • 

(8) 


These tensor elements correspond exactly to those of the NN RVB wavefunction. Note that the tensor elements of 
/ with the same indices identically vanish by symmetry. The solution of Eqs. (2) and (3) gives also non-zero tensor 
elements that correspond to quantum teleportation via diagonal bonds. 


R 


where Aq^ G R. One also gets (grouping the R and I tensors) 






ts _ 

^sss2 ^0 5 

(9) 

>s' _ \ss' 

^s2ss ^0 5 

(10) 

ts' _ \ss' 

^2sss ^0 5 

(11) 

}s' _ _ \ss' 

^ss2s ^0 ") 

(12) 

I tensors). 

— j- 4 

— A 2 -r 

(13) 

— _L 4 

— A 2 -T^A^hirab 

(14) 
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1/N 


1/N 


FIG. 5: Finite size scaling of the leading eigenvalues of the diagonal blocks (Fee and Foo) and of the off-diagonal blocks (fee 
and Feo) of the transfer matrix (after proper normalization). Same notation as Fig. [^b) in the main text. The three panels 
correspond to different choices of the parameters of the PEPS as indicated on the plot, (a) is a non-chiral RVB state (7ee = 1), 
(b) does not contain NN valence bonds and (c) is a repetition of Fig. [^b) for convenience. 


^s2ss 

^sss2 

^ss2s 

A^'- 

^2sss 

As'_ 

^ss2s 

^2sss 


■\ss 

^2 ~ ^^chirab 

(15) 

\ss' -xss 

^2 ^^chirab 

(16) 

^2 ^ ^^chirab 

(17) 

— A 2 -h ^A^hirab 

(18) 

ass' ;.\SS 

— A 2 — ^A^hirab 

(19) 

ass' n\ss' 

— ^2 ~ ^^chirab 

( 20 ) 


where G R are independent constants and s is the time-reversed of the spin s. 

Finally, enforcing the invariance under SU( 2 ) spin-rotations (singlet character of the PEPS) one finds that A^^ 
(^chiral (are) diagonal in the spin indices and even (odd) under spin inversion, 


Ar = ^..'A2, (21) 

Achiral ~ Achiral: (22) 

Ag^' = ^..^(-l)^Ao. (23) 

Note that spin rotation invariance leads also to the relation, 

Ao = 2 A 2 . (24) 


LEADING TM EIGENVALUES: COMPARISON BETWEEN DIFFERENT WAVE FUNCTIONS. 

We compare the finite size scalings of the LE of the TM for three choices of the tensor parameters in Fig. 
corresponding to the critical (non-chiral) RVB state in (a) and chiral topological liquids in (b) and (c). For all wave 
functions, the finite size scalings are consistent with 7^0 ^ 1 and 7 ^^^^ ^ 1. However, in (b) and (c) it is difficult to 
assess with full confidence that 77 ro/( 7 ee 77 r 7 r)^^^ and jeo/ ( 7 ee 7 oo)^^^ do not go to exactly 1 when Ny oc, as required 
to obtain separate topological sectors. 


ENTANGLEMENT SPECTRUM: COMPARISON BETWEEN DIFFERENT SIZES. 


Comparisons between ES for cylinders with different perimeters is shown in Fig. [^a-c). Remarkably the slope of 
the chiral edge modes is very similar for Ny = 6 and Ny = 8 ^_ 
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FIG. 6: Entanglement Spectrum vs momentum along the edge for Ny = 4,6,8 (for Ai = A 2 = Achiral = !)• (a), (b) and (c) 
show = 0, = 1 and IaS^I = 1/2, respectively. The dashed lines show a fit of the chiral edge modes for Ny = 8. (d) Finite 

size scaling of the lowest eigenvalue of the ES. 





















